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Abstract. The paper provides a review of A.M. Mathai's applications of the 
theory of special functions, particularly generalized hypergeometric functions, 
to problems in stellar physics and formation of structure in the Universe and 
to questions related to reaction, diffusion, and reaction-diffusion models. The 
essay also highlights Mathai's recent work on entropic, distributional, and dif- 
ferential pathways to basic concepts in statistical mechanics, making use of his 
earlier research results in information and statistical distribution theory. The 
results presented in the essay cover a period of time in Mathais research from 
1982 to 2008 and are all related to the thematic area of the gravitationally sta- 
bilized solar fusion reactor and fractional reaction-diffusion, taking into account 
concepts of non-extensive statistical mechanics. The time period referred to 
above coincides also with Mathai's exceptional contributions to the establish- 
ment and operation of the Centre for Mathematical Sciences, India, as well as the 
holding of the United Nations (UN)/European Space Agency (ESA)/National 
Aeronautics and Space Administration (NASA) of the United States/ Japanese 
Aerospace Exploration Agency (JAXA) Workshops on basic space science and 
the International Heliophysical Year 2007, around the world. Professor Mathai's 
contributions to the latter, since 1991, are a testimony for his social conscience 
applied to international scientific activity. 
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1 BIOGRAPHICAL NOTES: 
Professor Dr. A.M. Mathai 

Arakaparampil Mathai Mathai was born in Arakulam in the Idukki district of 
Kerala, India. He is the eldest son of Aley and Arakaparampil Mathai. He 
completed his High School education in 1953 with record marks in St. Thomas 
High School, Palai. He then entered St. Thomas College in Palai, obtained 
his BSc degree in mathematics in 1957 and went to the University of Kerala, 
Trivandruni, where he completed his MSc degree in statistics in 1959 with first 
class, first rank, and a golden medal. He joined St. Thomas College as a 
lecturer from 1959 to 1961. In 1961, he was awarded a Canadian Commonwealth 
Scholarship and wont to the University of Toronto wore he completed an MA 
and a PhD in mathematics and mathematical statistics in only three years, 
the minimum required time. He then joined McGill University in Montreal - 
one of the most prestigious universities in Canada - where he was Assistant 
Professor of Mathematics and Statistics, 1964-1968, Associate Professor, 1968- 
1978, and Pull Professor, 1979-2000. Since 2000, he is Emeritus Professor in 
Mathematics and Statistics at McGill University, Montreal, Canada. He served 
and is currently serving on editorial boards of a number of journals and he 
has supervised over 30 Master's students and 5 PhD students. Currently he 
is the Director of the Centre for Alatlicmatical Sciences, Trivandrum and Pala 
Campuses, Kerala, India which he is developing into an national, regional, and 
international centre of excellence for research in the mathematical sciences. The 
holding of Seven six-week SERC Schools in 1995, 2000, 2005, 2006, 2007, 2008, 
and 2009 are a testament of highest achievements in teaching and research in 
the mathematical sciences for the Centre. 

A.M. Mathai is a Fellow of the National Academy of Sciences, India, and a 
Fellow of the Institute of Mathematical Statistics, United States of America. 
He is an Elected Member of the International Statistical Institute, The Nether- 
lands. He is the Founder of the Canadian Journal of Statistics, Canada, and 
the Founder of the Statistical Science Association of Canada (now. Statistical 
Society of Canada). Professor Mathai is a Honorary Member of the Interna- 
tional Scientific Organizing Committee of the United Nations/European Space 
Agency /National Aeronautics and Space Administration of the United States of 
America/ Japanese Aerospace Exploration Agency Workshops on basic space sci- 
ence which have been held, since 1991, in India, Costa Rica, Colombia, Nigeria, 
Egypt, Sri Lanka, Germany, Honduras, Jordan, France, Mauritius, Argentina, 
China, United Arab Emirates, India, Japan, Bulgaria, and South Korea. Pro- 
fessor Mathai is author or co-author of fourteen research level books, seven 
undergraduate-graduate level books, and more than thirty edited and co-edited 
proceedings. More than 250 published research papers in journals from around 
the world are showing the breadth of the array of topics that he has tackled 
as well as the depth, the complexity and the originality of his published results 
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put him in a class all his own. All these publications concern A.M. Mathai's 
research areas with significant research contributions in mathematical statis- 
tics, applied statistics, probability, physics, mathematical physics, astrophysics, 
information theory, mathematics, and specific problems in biological models, 
orthogonal polynomials, dispersion theory, axiomatization of basic statistical 
concepts, pollution problems, and transportation problems. An electronic li- 
brary, recording all of A.M. Mathai's publications, is under preparation at the 
United Nations in Vienna. 

Around the world, A.M. Mathai is known for being very energetic, exceptionally 
industrious and remarkably gifted. Those individuals who have the good fortune 
to be close acquaintances of his are also aware of his moral integrity and his social 
conscience. The latter was particularly beneficial to have him as the advisor in 
the organization of scientific activities under the umbrella of the United Nations. 
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2 ASTROPHYSICAL THERMONUCLEAR 
FUNCTIONS 

Stars are gravitationally stabilized fusion reactors changing their chemical com- 
position while transforming light atomic nuclei into heavy ones. The atomic 
nuclei are supposed to be in thermal equilibrium with the ambient plasma 
(Maxwell-Boltzmann) . The majority of reactions among nuclei leading to a nu- 
clear transformation are inhibited by the necessity for the charged particles to 
tunnel through their mutual Coulomb barrier (Gamow). As theoretic;aI knowl- 
edge and experimental verification of nuclear cross sections increases it becomes 
possible to refine analytic representations for nuclear reaction rates. 

2.1 The Integrals 

The standard case of the thermonuclear function contains a nuclear cross sec- 
tion, the energy dependent term of a Taylor series, and the Maxwell-Boltzmann 
distribution function, 

hiz^u)"^ r y^e-ye-^y'Uy (1) 
Jo 

where y = E/kT and z — 27r(^/2fcT)^/^ZiZ2e^/?i. Considering dissipative 
collision processes in the thermonuclear plasma, cut off of the high energy tail 
of the Maxwell-Boltzmann distribution may occur, thus we write, 

l2{z,d,,.)=^[\^e-ye-^y~^dy, (2) 
Jo 

where d denotes a certain typically high energy threshold. 

Accommodating screening effects in the standard thermonuclear function 
we have to use nuclear cross section and the Maxwell-Boltzmann distribution 
function which leads to 

l3{z,t,u)=^ r y^e-ye-^^y+'^'Uy, (3) 
Jo 

where t is the electron screening parameter. 

Finally, if due to plasma effects a depletion of the Maxwell- Boltzmann 
distribution has to be taken into account, the thermonuclear function can be 
written in the following form 

h{z,5Ay)= y^e-ye-'y'e-^y'^dy, (4) 
Jo 

where the parameter S exhibits the enhancement or reduction of the high-energy 
tail of the Maxwell-Boltzmann distribution. 

Mathematical/statistical techniques can be used for deriving closed- form 
representations of the above thermonuclear functions. 
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2.2 Closed-form Representation of the Integrals 

For z > and i/ > 0, we have 

Jo o,i,i+.y 



For z > Q,d > Q, and > 0, we have 



^0 7r2 f-^ r! 



! \ M 



u+r+2 



i/+r+l,0,i. 
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We obtain 
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(11) 
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2.3 Further Reading 

Anderson, W.J., Haubold, H.J., and Mathai, A.M.: 1992, Astrophysical ther- 
monuclear functions, Astrophysics and Space Science, 214, 49-70. 

Bergstrom, Iguri, S., and Rubinstein, H.: 1999, Constraints on the variation 
of the fine structure constant from big bang nucleosynthesis. Physical Review 
D60, 045005. 

Critchfield, C.L.: 1972, Analytic forms of the thermonuclear function. In: Cos- 
mology, Fusion, and Other Matters. George Gamow Memorial Volume, Edited 
by F. Reines, University of Colorado Press, Colorado, pp. 186-191. 
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Vol. 348, Springer- Verlag, Berlin. 
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26, 400-413. 



3 SOLAR STRUCTURE IN TERMS OF 

GAUSS' HYPERGEOMETRIC FUNCTION 

Hydrostatic equilibrium and energy conservation determine the conditions in 
the gravitationally stabilized solar fusion reactor. We assume a matter density 
distribution varying non-linearly through the central region of the Sun. The 
analytic solutions of the differential equations of mass conservation, hydrostatic 
equilibrium, and energy conservation, together with the equation of state of 
the perfect gas and a nuclear energy generation rate e = eop"!"™, are given in 
terms of Gauss' hypergeometric function. This model for the structure of the 
Sun gives the run of density, mass, pressure, temperature, and nuclear energy 
generation through the central region of the Sun. Because of the assumption 
of a matter density distribution, the conditions of hydrostatic equilibrium and 
energy conservation are separated from the mode of energy transport in the 
Sun. 



3.1 Hydrostatic Equilibrium 

In the following we are concerned with the hydrostatic equilibrium of the purely 
gaseous spherical central region of the Sun generating energy by nuclear re- 
actions at a certain rate. For this gaseous sphere we assume that the matter 
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density varies non-linearly from the center outward, depending on two parame- 
ters 6 and 7, 

p{x) = pcfoix), (13) 

fn{x) = [1 - x'r, (14) 

where x denotes the dimensionless distance variable, x = r/i?0,O < a; < 1, Rq 
is the solar radius, 5 > 0, 7 > and 7 is kept a positive integer in the following 
considerations. The choice of the density distribution reveals immediately that 
p{x = 0) = pc is the central density of the configuration and p{x = 1) = is 
a boundary condition for hydrostatic equilibrium of the gaseous configuration. 
For the range < a; < 0.3 the density distribution can be fit numerically to 
computed data for solar models by choosing S = 1.28 and 7 = 10. The choice 
of restricting x to x < 0.3 is justified by looking at a Standard Solar Model 
which shows that a; < 0.3 comprises what is considered to be the gravitationally 
stabilized solar fusion reactor. More precisely, 95% of the solar luminosity is 
produced within the region x < 0.2(M < 0.3M©). The half-peak value for the 
matter density occurs at a; = 0.1 and the half-peak value for the temperature 
occurs at X = 0.25. The region x < 0.3 is also the place where the solar neu- 
trino fluxes are generated. As we are concerned with a spherically symmetrical 
distribution of matter, the mass M{x) within the radius x having the density 
distribution is 

M{x) = Mefuix), (15) 

where Mq denotes the solar mass and 2-P'i(-) is Gauss' hypergeometric function. 

This equation is satisfying the boundary condition M(x = 0) = and deter- 
mine the central value pc of the matter density through the boundary condition 
M{x = 1) = Mq, where pc depends then only on 6 and 7 of the chosen density 
distribution. 

For hydrostatic equilibrium of the gaseous configuration the internal pressure 
needs to balance the gravitational attraction. The pressure distribution follows 
by integration of the respective differential equation for hydrostatic equilibrium, 
making use of the above density distribution and mass distribution, that is 

q 



fpix) = 



(f + l)(f + 2)...(f+7) 



7! 



6^ 



E 



(-7). 



7! 



. 2Fi(-7, l+m;-+m+l;x^) 





,(18) 
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where G is Newton's constant and 2-Fi(-) denotes again Gauss' hypergeometric 
function. 

The Pochhammer symbol (| + m + 1)^ = r(| + m + 1 + 7)/r(| + m + 1) 
often appears in series expansions for hypergeometric functions. This equation 
gives the value of the pressure Pc at the centre of the gaseous configuration and 
satisfies the condition P{x = 1) = 0.. 



3.2 Equation of State 

It should be noted that P{x) denotes the total pressure of the gaseous configu- 
ration, that is the sum of the gas pressure and the radiation pressure (according 
to Stefan-Boltzmann's law) . However, the radiation pressure, although the ratio 
of radiation pressure to gas pressure increases towards the center of the Sun, 
remains negligibly small in comparison to the gas pressure. Thus, P{x) can be 
considered to represent the run of the gas pressure through the configuration 
under consideration. Further, the matter density is so low that at the temper- 
atures involved the material follows the equation of state of the perfect gas. 
Therefore, the temperature distribution throughout the gaseous configuration 
is given by 



r(| + l)(f + 2).-.(f +7) 



E 



(-7). 



[l-^'r^o^Kf+m)(|+m) 



(| + m+l)^ 



— X 



.6m+2 



3 



m: - +m + l\x ) 





(20) 



where k is the Boltzmann constant, Na Avogadro's number, /i the mean molec- 
ular weight, and 2-^1 (•) Gauss' hypergeometric function. The equations for 
T{x) reveal the central temperature for T{x = 0) = Tc and satisfy the bound- 
ary condition T{x = 1) = 0. Since the gas in the central region of the Sun 
can be treated as completely ionized, the mean molecular weight ^ is given by 
fj, = {2X + jY + ^Z)~^, where X,Y,Z are relative abundances by mass of 
hydrogen, helium, and heavy elements, respectively, and X + Y + Z = 1. 



3.3 Nuclecir Energy Generation Rate 

Hydrostatic equilibrium and energy conservation are determining the physical 
conditions in the central part of the Sun. In the preceeding Sections, the run of 
density, mass, pressure, and temperature have been given for a gaseous config- 
uration in hydrostatic equilibrium based on the equation of state of the perfect 
gas. In the following, a representation for the nuclear energy generation rate, 

e(p,T) = eop"T™, (21) 
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will be sought which takes into account the above given distributions of density 
and temperature and which can be used to integrate the differential equation of 
energy conservation throughout the gaseous configuration. In the above equa- 
tion, n denotes the density exponent, m the temperature exponent, and eo a 
positive constant determined by the specific reactions for the generation of nu- 
clear energy. 

3.4 Total Nuclecir Energy Generation 

The total net rate of nuclear energy generation is equal to the luminosity of the 
Sun, that means the generation of energy by nuclear reactions in the central 
part of the Sun has to continually replenish that energy radiated away at the 
surface. If L{x) denotes the outflow of energy across the spherical surface at 
distance x from the center, then in equilibrium the average energy production 
at distance x is 



Jo 

where p{x) and e{x) are given above. 
3.5 Further Reading 

Chandrasekhar, S.: 1939/1957, An Introduction to the Study of Stellar Struc- 
ture, Dover Publications Inc., New York. 

Clayton, D.D.: 1986, Solar structure without computers, American Journal of 

Physics, 54. 354-362. 

Haubold, H.J., Matliai, A.M.: 1995, Solar structure in terms of Gauss'hypergeometric 
function. Astrophysics and Space Science, 228, 77-86. 

Haubold, H.J.: 1995, An analytic solar model: Physical principles and math- 
ematical structures. International Journal of Mathematical and Statistical Sci- 
ences, 4, 31-41. 

Haubold, H.J., Mathai, A.M.: 2001, Sun, in Encyclopedia of Planetary Sciences, 
eds. J.H. Shirley and R.W. Fairbridge, Chapman and Hall, New York, pp. 786- 



Mathai, A.M. and Haubold, H.J.: 2008, Special Functions for Applied Physi- 
cists, Springer, New York. 

Mathai, A.M., Saxena, R.K., and Haubold, H.J.: 2009, The H-Function: Theory 
and Applications, Springer, New York. 

Stein, R.F.: 1966, in R.F. Stein and A.G.W. Cameron (eds.), Stellar Evolution, 
Plenum Press, New York. 




(22) 



794. 
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4 THE FRACTIONAL KINETIC EQUATION 



The subject of this remark is to derive the solution of the fractional kinetic 
equations. The results are obtained in a compact form containing the Mittag- 
LcfHcr fimction. which naturally occurs whenever one is dealing with fractional 
integral equations. 



4.1 Kinetic Equation and Mittag-Leffler Function 

In terms of Pochammer's symbol 

{a)n = {^'('; 

we can express the binomial series as 



,n=0 

x{a+l)...{a+n-l),neN 



a-)-" = E^- (23) 

The Mittag-LefHer function is defined by 

n— 

This function was defined and studied by Mittag-Leffler. We note that this 
function is a direct generalization of an exponential function, since 

Ei{z) :— exp{z). 

It also includes the error functions and other related functions, we have 

Ey2{±z^/^) = e^[l + erf{±z^l^)] = e'erfc{Tz^l% (25) 

where 

erf{z) ■=:^ erfc{z) := 1 - erf{z), z G C. (26) 

The equation 

gives a generalization of the Mittag-LcfFlcr fimction. When /3 = 1, the general- 
ization reduces to the original. Both functions are entire functions of order 1/a 
and type 1. The Laplace transform of Ea^/}{z) follows from the integral 

POO 

/ e-H^-^E„^f3{\ae)dt = p-^{l - ap-")-\ (28) 
Jo 
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where Re{p) > |a|i/", i?e(/3) > 0, which can be estabhshed by means of the 
Laplace integral 

/■oo 

/ e-PH''-'dt = Tip)/jf, (29) 
Jo 

where Re{p) > 0, Re{p) > 0. The Riemann-Liouville operator of fractional 
integration is defined as 



aD-'^m = f^J^ fi^Kt - uy-^du, v>0, (30) 
with aD^f{t) = f{t). By integrating the standard kinetic equation 

j^Ni{t) = -CiNi{t),{ci>0), (31) 

it is derived that 

N,it) -No = -c, oD-^N,it), (32) 

where o^T^ the standard Ricniann integral operator. The number density of 
species i,Ni = Ni{t), is a function of time and Ni{t = 0) = Nq is the number 
density of species i at time t = 0. By dropping the index i, the solution of its 
generalized form 

N{t) -No = -c" oD^'Nit), (33) 

is obtained as 

By virtue of the original definition of the Mittag-Lcfflcr function wc can rewrite 
this solution in terms of the Mittag-Lcfflcr function in a compact form as 

N{t) = NoE^i-c^t"), u>Q. (35) 



4.2 Conclusions 

Solutions are given in terms of the ordinary Mittag-LefHer function and their 
generalization, which can also be represented as FOX's H-functions. The or- 
dinary and generalized Mittag-Leffler functions interpolate between a purely 
exponential law and power-like behavior of phenomena governed by ordinary 
kinetic equations and their fractional counterparts. 



4.3 Further Reading 
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diffusion equations. Astrophysics and Space Science, 305, 289-296. 
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5 GENERALIZED MEASURE OF ENTROPY, 
MATHAI PATHWAYS, AND TSALLIS 
STATISTICS 

When deriving or fitting models for data from physical experiments very often 

the practice is to take a member from a parametric family of distributions. But 
it is often found that the model requires a distribution with a more specific 
tail than the ones available from the parametric family, or a situation of right 
tail cut-off. The model reveals that the underlying distribution is in between 
two parametric families of distributions. In order to create a distributional 
pathway for proceeding from one functional form to another a pathway param- 
eter is introduced and a pathway model is created. This model enables one to 
go from a generalized type-1 beta model to a generalized type-2 beta model 
to a generalized gamma model when the variable is restricted to be positive. 
More families arc available when the variable is allowed to vary over the real 
line. Corresponding multitudes of families are available when the variable is in 
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the complex domain. Originally, Mathai deals mainly with rectangular matrix- 
variate distributions in the real case and the scalar case is the topic of discussion 
in the current remark because the scalar case is connected to many problems in 
physics. Similarly, through a generalized measure of entropy, an entropic path- 
way can be devised that leads to different families of entropic functionals that 
produce families of distributions by applying the maximum entropy principle. 



5.1 The Pathway Model in the Real SceJar Case 

For the real scalar case the pathway model for positive random variables is the 
following: 

f{x) = cx^^^[l-a{l-a)x^]T^, x > 0, (36) 

a>0,(5>0,/3>0,l — a(l — a)x^ > 0, 7 > where c is the normalizing constant 
and a is the pathway parameter. The normalizing constant in this case is the 
following: 



r(i)r(T^ + 
r(i)r(/T-i) 

S (a/3)T 



for a < 1 (37) 



for - ^ > 0, a > 1 (38) 

a — 1 



for a^l. (39) 



5.2 Pathway Model from a Generalized Entropy Measure 

A generalized entropy measure of order a is a generalization of Shannon en- 
tropy and it is a variant of the generalized entropy of order a. In the dis- 
crete case the measure is the following: Consider a multinomial population 
P = (Pi) • • • )Pfe)) Pi > 0, z = 1, . . . , fc, pi -h . . . +Pfc = 1. Define the function 

MkAP) = " ~ \ Q 7^ 1, -00 < a < 2 (40) 

a — 1 

k 

lim Mfe,„(P) = - Vft Inpi = Sk{P) (41) 

a— >1 

j=l 

by using L'Hospital's rule. In this notation In is taken as zero when any 
Pi = 0. Thus it is a generalization of Shannon entropy Sk{P)- Note that it is a 
variant of Havrda-Charvat entropy Hk,a{P) and Tsallis entropy Tk^a{P) where 

^kAP) = ^|^-'cf[7^ a 7^ 1, a > (42) 
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and ^ 

TkAP)= ^'f^'~'^ , ">0. (43) 

1 — a 

We will introduce another measure associated with the above and parallel to 
Renyi entropy Rk,a in the following form: 



2-a 
= 1 f'i 



In (Ell p: 

Renyi entropy is given by 



Ml„{P) = ^ a^l,-oo<a<2. (44) 



In 



(E-=iPr) 



RkAP)= \ a^l,a>Q. (45) 



5.2.1 Continuous analogue 

The continuous analogue to the generalized entropy measure is the following: 

M„(/) = ^-'^''^ ^ (46) 

a — 1 

= ^1 = ^1 '"^l'"<2, 

where E[-] denotes the expected value of [•]. Note that when a = 1, 
E[f{x)Y- = E[f{x)]° = 1. 

It is easy to see that the generalized entropy measure is connected to Kerridge's 
"inaccuracy" measure. The generalized inaccuracy measure is E[q{x)Y~'^ where 
the experimenter has assigned q{x) for the true density f{x), where q{x) could 
be an estimate of f{x) or q{x) could be coming from observations. Through 
disturbance or distortion if the experimenter assigns [f{x)Y~'^ for [^(a;)]"'^"" 
then the inaccuracy measure is Ma{f). 



5.2.2 Distributions with maximum generalized entropy 

Among all densities, which one will give a maximum value for Ma(/)? Con- 
sider all possible functions f{x) such that f{x) > for all x, f{x) = outside 
(a, h), a < b, f{a) is the same for all such f{x), f{b) is the same for all such /, 
f{x)dx < 00. Let f{x) be a continuous function of x with continuous deriva- 
tives in (a, b). Let us maximize jl^[f{x)]'^^"dx for fixed a and over all functional 
/, under the conditions that the following two moment-like expressions be fixed 
quantities: 

I x(''-i)(i-")/(a;)da; = given, and / x^''-^'^^^-'''^+^ f{x)dx = given (47) 

J a J a 
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for fixed 7 > and d > 0. Consider 

U = [/(a;)]2-« - Aia;('>'-i)(i-")/(a;) + X2X^'''-'^^^^-"^+^ f{x), a < 2, a ^ 1 (48) 

wlicrc Ai and A2 are Lagrangian multipliers. Then the Euler equation is the 
following: 

^=0 ^ (2-a)[/(a;)]i-«-Aia;(T-i)(i-«)+A2ar(T-i)(i-«)+* = 

^ [/(^)]^-" = 7^-(^-^)^^-"^[l-^^'] (49) 
(2 — a) Ai 

^ /(a;) = c a;T-i [1 - r/(l - a)x''] ^ (50) 

where A1/A2 is written as 77(1 — a) with 77 > such that 1 — 77(1 — a)x^ > 
since f{x) is assumed to be non-negative. By using the conditions above we can 
determine c and 77. When the range of x for which f{x) is nonzero is (0, 00) and 
when c is a normalizing constant, then f{x) is the pathway model of Mathai in 
the scalar case where a is the pathway parameter. When 7 = 1,(5 = 1 in f{x) 
then f{x) produces the power law. The form above for various values of Ai and 
A2 can produce all the four forms 

aia;^-^[l - /3i(l - a)x^]-T^, a2a;^"^[l - /?2(1 - a)a;*]T^fora < 1 

and 

aax'^-^ll + (33{a - l)a;'']"^, a^x^^-^ll + f3i{a - l)x'']^fora > 1 

with aijPi > 0,i = 1,2,3,4. But out of these, the second and the third forms 
can produce densities in (0, 00). The first and fourth will not be converging. 
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6 GRAVITATIONAL INSTABILITY IN A 
MULTICOMPONENT COSMOLOGICAL 
MEDIUM 

The paper presents results for deriving closed-form analytic solutions of the 
non-relativistic linear perturbation equations, which govern the evolution of 
inhomogeneities in a homogeneous spatially flat multicomponent cosmological 
model. Mathematical methods to derive computable forms of the perturbations 
are outlined. 

6.1 Introduction 

These closed-form solutions are valid for irregularities on scales smaller than 

the horizon. They may be used for analytical interpolation between known ex- 
pressions for short and long wavelength perturbations. We present the physical 
parameters relevant to the problem and the fundamental differential equations 
for arbitrary polytropic index "fi and a radiation or matter dominated Universe, 
respectively. Solutions of this fundamental equation are provided for relevant 
parameter sets including the polytropic index 7j and the expansion law index rj. 
To catalogue these solutions in closed-form, the theory of Meijer's G- function 
will be employed. 

6.2 Some Parameters of Physical Significance 

The growth of density inhomogeneities can begin as soon as the Universe is 
matter-dominated. Baryonic inhomogeneities cannot begin to grow until after 
decoupling because until then, baryons are tightly coupled to the photons. After 
decoupling, when the Universe is matter-dominated and baryons are free of the 
pressure support provided by photons, density inhomogeneities in the baryons 
and any other matter components can grow. Actually the time of matter- 
radiation equality is the initial epoch for structure formation. In order to fill 
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in the details of structure formation one needs "initial data" for that epoch. 
The initial data required include, among others, (1) the total amount of non- 
relativistic matter in the Universe, quantified by Qq, and (2) the composition of 
the Universe, as quantified by the fraction of critical density, fli = Pi/pc, con- 
tributed by various components of primordial density perturbations (i= baryons. 
Wimps, rclativistic particles, etc.). Here the critical density pc is the total mat- 
ter density of the Einstein-de Sitter Universe. Speculations about the earliest 
history of the Universe have provided hints as to the appropriate initial data: 
fio = 1 from inflation; 0.015 < fis < 0.15 and ^Iwimp ^ 0.9 from inflation, 
primordial nucleosynthesis, and dynamical arguments. In the following we as- 
sume that X^Oj = l,Oi = const.. The cosmological medium is considered 
to be unbounded and that there may be a uniform background of rclativistic 
matter. It is further assumed that the matter is only slightly perturbed from 
the background cosmological model. This assumption may give a good descrip- 
tion of the behavior of matter on large scales even when there may be strongly 
nonlinear clustering on small scales (primordial objects). Also assumed is that 
matter can be approximated as an ideal fluid with pressure a function of density 
alone. It consists of i components having the densities pi and velocities of sound 
Pf = dPi/dpi (X. p~'^~^ , when an equation of state 

Pi oc p^- (51) 

has been taken into account. The i components of the medium arc interrelated 
through Newton's field equation = ^-nGY^pi, containing the combined den- 
sity ^ pi of all components. Superimposed upon an initially homogeneous and 
stationary mass distribution shall be a small perturbation, represented by a sum 
of plane waves 

5=^-^ = 5i{t)e'''\ (52) 
P 

where k = 27ra/A defines the wave number k and A is the proper wavelength. In 
the linear approach the system of second order differential equations describing 
the evolution of the perturbation in the non-relativistic component i is 

+ 2(^)^ + k'^m = 47rG^ft5„ i = 1, ...,m, (53) 

i=i 

where for all perturbations the same wave number k is used. This equation is 
valid for all sub-horizon-sized perturbations in any non-relativistic species, so 
long as the usual Friedmann equation for the expansion rate is used. Although 
the following investigation of closed-form solutions of the fundamental differen- 
tial equation governing the evolution of inhomogeneities in a multicomponent 
cosmological medium is quite general, it will be presented within the context of 
the inflationary scenario. Therefore in the following we consider an Einstein-de 
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Sitter Universe with zero cosmological constant. After inflation the Friedmann 
equation for the cosmological evolution reduces to 

i-r = - —p. (54) 

The continuity equation gives p oc and one has a oc i^l^ and pi = Q.i/{QiTGt^). 
During the expansion the Hubble parameter changes as H = r]t~^, where the 
expansion law index is 77 = ^ in the radiation-dominated epoch and 77 = | in 
the matter-dominated epoch respectively. The wave number is proportional to 
so that = kft'^^^~'^~'^^\ where now the constants ki come from both 
the wave vector and the velocity of sound. If ki = the sound velocity f3i equals 
zero, in which case the adiabatic index ji loses its sense. Defining the param- 
eter a, = 2(2 — T] — ji) that absorbs the adiabatic index as well as the type of 
expansion law we can write: 

2 

t^6i{t) + 2r]t6^{t) + kle^Si = 3 XI ^J'^j(^)' i=l,---,m, (55) 

and dots denote derivatives with respect to time. We introduce a time oper- 
ator A = M/dt and change the dependent variable, the density perturbation 
5i(t), dealing with the function <l>j instead of 5i by setting 6i{i) — t"^i{t). The 
equation for the function $j is then given by 

2 " 



where 



$i = t-"6i, h = kft^^ - a^, a, = 2(2 - 77 - 7,), a = -[ ) . (56) 



Observing that T,^j = 1 and operating on both sides of the differential equation 
for the $j by we obtain the fundamental equation 

c) n rn 

I^''^i + A\hi^i)--{A^^i + hi^i) = --Y,hjnj^j, i = l,...m. (57) 

i=i 

As indicated above for 77 the values | and ^ are significant and in what follows 
it will be assumed that 2 > 74 > | . We have chosen the range of values of 7j for 
physical as well as mathematical reasons as will be evident later in the analysis. 
For some of these parameter values we will consider a multicomponent medium. 
Consider the special case 

6r = 

62 = 63 = ... = 6„ = 6 = fcV--^^^^^^ (58) 
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of the fundamental equation. Let k = 0. For this case the parameters are the 
following: 



ai,a2 = -1 ± 



,bl,bl = ± 
(2,7-1)2 



(2^-^ 2 
4 3 



2 2^ 



(59) 



6.3 Solution when the Parameters Differ by Integers 

The general solution are for finite values of t and for the cases that the b*'s do 

3) 



one has b^ — 63 



-1 and bl — b^ 



-1. 



not differ by integers. At two points, namely (77 = |, 7j = |) and (77 = |, 7j = 1) 
6.4 Modifications of the G'jS for the Cases (77 = i)7i — 1) 



In 



and {rj = 

Consider, for x = kft"'^ / al, ai ^ 0, 

r(i + s)r(-al - s)T{-ai - s)x-^ds 



r(t-s)r(-i-s)r(|-s) 



(60) 



Note that a zero coming from r(— ^ — s) at s = — | coincides with a pole coming 
from r{j + s) at s = — J. This can be removed by rewriting as follows: 

r(| + s)r{-al - s)r{-a*2 - s)x-^ds 



Gi = 



1 

"2^1 



r(|-s)r(f-5)r(|-s) 



(61) 



Evaluating this expression as the sum of the residues at the poles of r(| + s) 

,r(-at + |)r(-a^ + |) 



one has 



r(f)r(2)r(M) 

r. / * 5 , 5 10 ^ 14 
X 2i^3(-ai + ^, + ^; ^, 2, — ; -x). 



(62) 



Evaluating G4 also the same way one has the following: 

r{-^ + s)T{-at-s)T{-a*2-s) 



-I 
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r r{-l + s)r{-al-s)T{-a*^-s) 

^TriJL r(|-s)r(|-s)r(-i-s) 

T{-k-al)r{-\-a*,) 

X 4 



X 



ds 



r(i)r(2)r(-i) 




(63) 



Thus the complete solution in this case is of the form 



$1 = Cid + C2G2 + C3G3 + C4G4 



(64) 



where ci, C2, C3, a are arbitrary constants. 

6.5 Conclusion 

We have presented closed-form solutions of the non-rclativistic linear perturba- 
tion equations which govern the evolution of inhomogeneities in a spatially flat 
multicomponent cosmological medium. The general solutions are catalogued 
according to the polytropic index 7^ and the expansion law index r] of the mul- 
ticomponent medium . All general solutions are expressed in terms of Meijer's 
G-function and their Mellin-Barnes integral representation. The proper use of 
this function simplifies the derivation of solutions of the linear perturbation 
equations and opens ways for its numerical computation. 
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7 Developing Basic Space Science World-Wide 
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the workshops will develop the International Space Weather Initiative. 

7.1 Further Reading 

Wamstcker, W., Albrccht, R., and Haubold, H.J. (Eds.): 2004, Developing 
Space Science World-Wide: A Decade of UN/ES A Workshops, Kluwer Aca- 
demic Publishers, Dordrecht Boston London. 

Thompson, B.J., Gopalswamy, N., Davila, J.M., and Haubold, H.J. (Eds.): 

2009, Putting the "F in IHY: The United Nations Report for the International 
Heliophysical Year 2007, Springer, New York. 

Briand, C, Antonucci, E., Gopalswamy, N., Eichhorn, C, Sakurai, T., Mathai, 
A.M., and Haubold, H.J. (Eds.): 2009, International Heliophysical Year 2007, 
Second European General Assembly, Italy, and Third UN/ESA/NASA Work- 
shop, Japan. Earth, Moon, and Planets, VoL 104, 1-360. 



21 



